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Abstract 

The perturbative QCD predicts that the growth of the gluon density at small- x (high energies) 
should saturate, forming a Color Glass Condensate (CGC), which is described in mean field approx- 
imation by the Balitsky-Kovchegov (BK) equation. Recently, the next-to-leading order corrections 
for the BK equation were derived and a global fit of the inclusive ep HERA data was performed, 
resulting in a parameterization for the forward scattering amplitude. In this paper we compare 
this parameterization with the predictions of other phenomenological models and investigate the 
saturation physics in diffractive deep inelastic electron-proton scattering and in the forward hadron 
production in pp collisions. Our results demonstrate that the running coupling BK solution is able 
to describe these observables. 
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I. INTRODUCTION 



The understanding of the high energy (small x) regime of Quantum Chromodynamics 
(QCD) has been one of the main challenges of this theory, which has been intensely in- 
vestigated through high energy collision experiments. This regime, where one expects to 
observe the non-linear behavior predicted by theoretical developments, has been explored in 
ep collisions at DESY-HERA and pp/dA collisions at BNL-RHIC and, in a near future, in 
pp/pA/AA collisions at CERN-LHC. In particular, at high energies, the growth of the par- 
ton distribution is expected to saturate, forming a Color Glass Condensate (CGC), whose 
evolution with energy is described by an infinite hierarchy of coupled equations for the 
correlators of Wilson lines 

q y, s, q, q, y, 0, y, y, y y y y . 

In the mean field approx- 



imation, the first equation of this hierarchy decouples and boils down to a single non-linear 
integro-differential equation: the Balitsky-Kovchegov (BK) equation 
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,H5[. 



The BK equation determines, in the large- N c (the number of colors) limit, the evolution of 
the two-point correlation function, which corresponds to the scattering amplitude M(x, r, b) 
of a dipole off the CGC, where r is the dipole size and b the impact parameter. This quantity 



encodes the information about the hadronic scattering and then about the non- 
quantum effects in the hadron wave function (For recent reviews, see e.g. 

y 
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inear and 
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id). 



In the last years, several groups have constructed phenomenological models which satisfy 



the asymptotic be 



RHIC data 120 



raviors of the leading order B 
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equation in order to fit the HERA and 
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361 ] . In general, it is 



assumed that the impact parameter dependence of M can be factorized as M^^{x,r,b) = 
J\fjr A (x,r)S(b), where S(b) is the profile function in impact parameter space and M? and 
Ma are the fundamental and adjoint dipole scattering amplitudes, respectively. The latter 
can be modelled in the coordinate space, through a simple Glauber-like formula, which reads 

1 



M A (x,r) 



cxp 



(i) 



where 7 is the anomalous dimension of the target gluon distribution, and the former can be 
also parameterized as in ([T]), after the replacement Q 2 S — > Q 2 s Cf/Ca = 4/9 Q^. The main 
difference among the distinct phenomenological models comes from the behavior predicted 
for the anomalous dimension, which determines the transition from the non-linear to the 
extended geometric scaling regime, as well as from the extended geometric scaling to the 
DGLAP regime. In this paper we restrict our analyses to the model proposed in Ref. [36j , the 



so called BUW model, which is able to describe the ep HERA data for the proton structure 



IC energy 
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34j models, which 



function and the hadron spectra measured in pp and dAu collisions at RE 
It should be contrasted with, for instance, the IIM jsjj] and DEIJ 
are only able to describe one of these data sets. Another feature of the BUW model which 
motivates this analysis is that it explicitly satisfies the property of geometric scaling, which 
is predicted for the solutions of the BK equation in the asymptotic regime of large energies. 
In the BUW model, the adjoint dipole scattering amplitude is parameterized in momentum 
space and is given by 

1 



J\f A {x,p T ) 



d re 



(2) 



where 7 is assumed to be a function of pt, rather than r, in order to make easier the 
evaluation of its Fourier transform, and is given by 7(j»t,x) = 7 S + A^y{pT,x), where 7 S = 



0.628 and 



36] 



A7(j9 T , x) = Ajbuw = (1 - 7 S 



(3) 



{u a - 1) + b' 

In the expression above, u = Pt/Qs{%) an d the two free parameters a = 2.82 and b = 168 



are fitted in such a way to describe the RHIC data on hadron production 
Eq.flSD, that this model satisfies the property of geometric scaling [37], 
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t is c 
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ear, from 



401 ]. since 



A7 depends on x and pr only through the variable pt/Q s {x). Besides, in comparison with 
other phenomeno logical parameterizations, in the BUW model, the large pt limit, 7 — > 1, 



is approached much faster, which implies different predictions 



hadron and photon yield (For a detailed discussion see Ref. 411]) . 



or the large pr slope of the 
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Recently the next-to- leading order corrections to BK equation were calculated 42j, |43, 



45 



471 ] through the ressumation of a s Nf contributions to all orders, where Nf is 



the number of flavors. Such calculation allows one to estimate the soft gluon emission 
and running coupling corrections to the evolution kernel and, in particular, the authors 
have verified that the dominant contributions come from the running coupling corrections, 



which allows to determine the scale of the running coupling in the 
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ternel. The solution 
Basically, one has 



of the improved BK equation was studied in detail in Refs. 
that the running of the coupling reduces the speed of the evolution to values compatible 
with experimental data, with the geometric scaling regime being reached only at ultra-high 



energies. In Ref. 



481 ] . the solution of the improved BK equation was used to calculate 



the pseudorapidity density of charged particles produced in nucleus-nucleus collisions and a 
remarkable good agreement with the RHIC data was observed. 

More recently, a global analysis of the small x data for the proton structure function 



using the improved BK equation was performed |49j (See also Ref. 55]). In contrast to the 
BK equation at leading logarithmic a s \n(l/x) approximation, which fails to describe data, 
the inclusion of running coupling effects to evolution renders BK equation compatible with 
them. The impact parameter dependence was not taken into account, the normalization of 
the dipole cross section was fitted to data and two distinct initial conditions, inspired in the 
Golec Biernat-Wusthoff (GBW) 20j and McLerran-Venugopalan (MV) 56j models, were 
considered. The predictions resulted to be almost independent of the initial conditions and, 
besides, it was observed that it is impossible to describe the experimental data using only 
the linear limit of the BK equation, which is equivalent to Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) equation [57 1. 



In this paper we compare the parameterization proposed in |49|] with the predictions of 
the BUW model. A basic difference among these models is the set of experimental data used 
in order to constrain the free parameters of the parameterizations or the initial condition 
of the running coupling (RC) BK evolution equation. For instance, while the BUW model 



361 ] uses the RHIC data for the forward hadron production in dAu collisions, in Ref. 49] 



the initial condition is constrained using the small-x HERA data for the proton structure 
function. In the original papers, the predictions of these two models were compared with the 
cr 7 * p and Fl HERA data, respectively. However, the resulting predictions for observables 
measured in diffractive deep inelastic scattering (DDIS) are still an open question. This is 
one of our goals. Another one is to calculate for the first time the forward hadron production 
in pp collisions using the solution of the RC BK evolution equation. This allows one to check 
the generalized dilute-dense factorization used to calculate the observables in hadron-hadron 
collisions and to verify if it is possible an unified description of the high energy regime in ep 
and pp collisions using the Color Glass Condensate formalism. 

This paper is organized as follows. In Section II we present the BK equation in both 
leading-order and with the improvements provided by the inclusion of the running of the 
coupling, and compare the behavior of the resulting dipole-target scattering amplitude with 
well-known parametrizations based on the saturation physics. Section III is devoted to the 
description of recent HERA data on the diffractive proton structure function F^ 3 \ while 
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in Section IV hadron production in pp collisions at RHIC is investigated. The conclusions 
are presented in Section V. 



II. BK EQUATION AND RUNNING COUPLING EFFECTS 

The Balitsky-Kovchegov equation is the simplest non-linear evolution equation for the 
dipole-hadron scattering amplitude, being actually a mean field version of the first equation 
of the Balitsky hierarchy . In leading order (LO), and in the translational invariance 
approximation — in which the scattering amplitude does not depend on the collision impact 
parameter b — it reads 



dX F (r,Y) 



dY 



J dr ± K h °(r, r 1? r 2 )[M F (n, Y)+Af F (r 2 , Y)-Af F (r, Y)-Af F ( n , Y)M F (r 2 , Y)], 

(4) 

where M F {r,Y) is the (fundamental) scattering amplitude for a dipole (a quark- ant iquark 
pair) off a target, with transverse size r = \r\, Y = \r(xq/x) (xq is the value of x where the 



evolution starts), and r 2 = r — T\. K LO is the evolution kernel, given by 

K LU r,r 1? r 2 = -— , 5 

2tt r(r 2 

where a s is the (fixed) strong coupling constant. This equation is a generalization of the 
linear BFKL equation (which corresponds of the first three terms), with the inclusion of the 
(non-linear) quadratic term, which damps the indefinite growth of the amplitude with energy 

re same universality class of 



591 ] and, as a consequence, it 



predicted by BFKL evolution. It has been shown [58j to be in t 
the Fisher-Kolmogorov-Pertovsky-Piscounov (FKPP) equation 
admits the so-called traveling wave solutions. This means that, at asymptotic rapidities, the 
scattering amplitude is a wavefront which travels to larger values of r as Y increases, keeping 
its shape unchanged. Thus, in such asymptotic regime, instead of depending separately on 
r and Y, the amplitude depends on the combined variable rQ s (Y), where Q S (Y) is the 
saturation scale. This property of the solution of BK equation is a natural explanation to 
the geometric scaling, a phenomenological feature observed at the DESY ep collider HERA, 



in the measurements of inclusive and exclusive processes 33, |38|, |39|, |40 ] . 

Although having its properties been intensely studied and understood, both numerically 
and analytically, the LO BK equation presents some difficulties when applied to study DIS 



small-x data. In particular, some studies concerning this equation |50|, |51 



152, 153, 154J have 

shown that the resulting saturation scale grows much faster with increasing energy (Q 2 S ~ 



x 



with A ~ 4.88A c a s /7r rs 0.5 for a s = 0.2) than that extracted from phenomenology 



(A ~ 0.2 — 0.3). This difficulty could be solved by considering smaller values of the strong 
coupling constant a s , but this procedure would lead to physically unrealistic values. One 
can conclude that higher order corrections to LO BK equation should be taken into account 
to make it able to describe the available small- a; data. 

The calculation of the running coupling corrections to BK evolution kernel was explicitly 



performed in 42 



43 



45 



46 



471 ]. where the authors included a s Nf corrections to the kernel 



to all orders. The improved BK equation is given in terms of a running coupling and a 
subtraction term, with the latter accounting for conformal, non running coupling contribu- 
tions. In the prescription proposed by Balitsky in [45] to single out the ultra-violet divergent 
contributions from the finite ones that originate after the resummation o f q uark loops, the 
contribution of the subtraction term is minimized at large energies. In [49] this contribu- 
tion was disregarded, and the improved BK equation was numerically solved replacing the 
leading order kernel in Eq. (jU) by the modified kernel which includes the running coupling 



corrections and is given by [4| 

N c a s (r 2 ) 



K B *\r,r u r 2 ) 



2tt 2 



r 2 If a s {r\) 



1 +^ 



1 f a s (r 2 ) 



- 1 



(6) 



44j, it has been confirmed that 



_r\r\ r\ \a s (r%) 

From a recent numerical study of the improved BK equation 
the running coupling corrections lead to a considerable increase in the anomalous dimension 
and to a slow-down of the evolution speed, which implies, for example, a slower growth of 
the saturation scale with energy, in contrast with the faster growth predicted by the LO 
BK equation. As it was discussed in the Introduction, the improved BK equation has been 
shown to be really successful when applied to the description of the ep HERA data for the 
proton structure function, which motivates us to extend the study for other observables. 



A. BK with running coupling and phenomenological models 

Before performing our phenomenological study using the improved BK equation, whose 
solution we will call RC BK from now on, it is interesting, at this point, to investigate the 
behavior of its solution and compare it with those from the phenomenological models based 
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FIG. 1: Dependence of the fundamental dipole scattering amplitude in the squared pair separation 
r 2 at different values of x. 

on saturation physics, which are able to describe HERA and/or RHIC data. It is important 
to point out that in order to make this comparison, as well as the phenomenological study 
;o be developed in the following sections, we make use of the public- use code available in 



60]. 



In Fig{T]we show the pair separation dependence of the fundamental scattering amplitude 
A/jr for different values of x (For a related discussion see Ref. [61]). From this figure one can 



32] and BUW parameterizations present a similar behavior 



observe that while the KKT 
for small r 2 , the RC BK one predicts a smoother dependence. In the limit of large pair 



separations, or large dipoles, the IIM 31| parameterization saturates, while the BUW, KKT 
and RC BK ones still present a residual dependence, demonstrating that for these models, 
the asymptotic regime is only reached for very large pair separations. A characteristic 
feature which is evident in the IIM model is that the dipole scattering amplitude saturates 
for smaller dipoles when x assumes smaller values. An important aspect to be emphasized 
is the large difference between the predictions in the transition region, which is expected to 
be probed at HERA and RHIC. 

A basic feature of the phenomenological parameterizations is that they grow in the region 
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FIG. 2: Energy dependence of the fundamental scattering amplitude as a function of x at different 
photon virtualities. 

of small values of rQ s as a power of rQ s , i.e. Af(r, x) oc {rQ s ) 2 ~ <£ ff . However, j e ff is different 
in each model, being 7 < 1 for the IIM model and about | for the KKT and BUW ones. 
This implies a different rQ s dependence of the dipole scattering amplitudes and dipole cross 
sections. Since the saturation scale drives the energy dependence of the dipole cross section, 
these models present a very distinct energy dependence. This can clearly be seen in Fig. [5J, 
where we present the x dependence of the dipole scattering amplitudes for different values 
of the squared pair separation given by r 2 = 1/Q 2 . We observe that for large Q 2 (small 
pair separation) the dipole scattering amplitude is dominated by the linear limit. Since the 
models have different behaviors in this limit, the energy dependence is also different, with 
the IIM model presenting the strongest growth at small x. On the other hand, the KKT and 
BUW models predict the two smallest growths with the energy. At large pair separations 
r > 1/Qs, which characterize the saturation regime, the IIM model predicts the saturation 
of the dipole scattering amplitude, while the KKT and BUW ones still present a growth 
at small values x. In comparison, the RC BK solution predicts a smooth growth, with a 
smaller normalization and a delayed saturation of the fundamental scattering amplitude at 
small values of x. Basically, the asymptotic saturation regime is only observed for very small 
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values of x, beyond the kinematical range of HERA. 



III. DIFFRACTIVE DIS IN THE CGC FORMALISM 



In the last years, significant progress in understanding diffraction has been made at the 
ep collider HERA (See, e.g. Refs. [62|, |63|, |64(). Currently, there exist many attempts to 



describe the diffractive part of the deep inelastic cross section within pQCD (See, e. g. R efs. 



20, 



65, 



66 



67j). One of the most successful approaches is the saturation one [2l|, |29j, [65] 



based on the dipole picture of DIS 



68 



691 ] . It naturally incorporates the description of 



both inclusive and diffractive events in a common theoretical framework, as the same dipole 
scattering amplitude enters in the formulation of the inclusive and diffractive cross sections. 
In the studies of saturation effects in DDIS, non- linear evolution equations for the dipole 

72], new measurements proposed 0, Q, 



scattering amplitude have been derived 70, 

HQ 



71 



77 



781 ] . As shown in Ref. 



21], the total 



and the charm contribution estimated 
diffractive cross section is much more sensitive to large-size dipoles than the inclusive one. 
Saturation effects screen large-size dipole (soft) contributions, so that a fairly large fraction 
of the cross section is hard and hence eligible for a perturbative treatment. Therefore, the 
study of diffractive processes becomes fundamental in order to constrain the QCD dynamics 
at high energies. 



The diffractive process can be analyzed in detai 



structure function F2 (Q 2 , (3, xjp) . In Refs. [21 



68 



studying the behavior of the diffractive 
] the authors have derived expressions for 



,D(3) 



directly in the transverse momentum space and then transformed to impact parameter 



space where the dipole approach can be applied. Following Ref. 2l| we assume that the 
diffractive structure function is given by 



-,D(3) 



(Q\(3 } xjp) = F q l L + F n D , T + F D 



qq,T 



qqg? 



(7) 



where T and L refer to the polarization of the virtual photon. For the qqg contribution 
only the transverse polarization is considered, since the longitudinal counterpart has no 
eading logarithm in Q 2 . The computation of the different contributions was done in Refs. 



21 



68 



79] and here we quote only the final results: 

3Q 6 



*l/2 



xjpF»(Q\i3,xjp) 



32n 4 (3B D 



r 

J^ej2 / daa\l - a) 3 $ , 



(8) 



x F F q l T (Q 2 ,P,xjp) 



3Q 4 



l28n 4 (3B D 



2J ej2 / daa(l - a) {e 2 [a 2 + (1 - a) 2 ]$i + mj$ } 



(9) 



where the lower limit of the integral over a is given by a>o 



1 — jtff I wo 



have introduced the auxiliary functions 651 ] : 

poo 

$o,i=( / rdrK ^(er)a d i P (x]p,r)J ^(kr] 



(10) 



Following Refs. 



21 



78 



79 



801 ] . here we calculate the qqg contribution within the dipole 



picture at leading In Q 2 accuracy, where it reads 



dk 2 \n 



iiPats 



dz 



i-2 

z 



k 2 t 



l-zf 
udu a dip (u/k t ,x I p)K2 



(11) 



-ir 



JM 



As pointed in Ref. 



291 ]. at small (3 and low Q 2 , the leading m(l//3) terms should be resumed 



and the above expression should be modified. However, as a description with the same 
quality using the Eq. ffTTl) is possible by adjusting the coupling j^, in what follows we will 
use this expression for our phenomenological studies. We use the standard notation for the 
variables (3 = Q 2 /(M 2 X + Q 2 ), x F = (M| + Q 2 )/{W 2 + Q 2 ) and x = Q 2 /{W 2 + Q 2 ) = pxjp, 
where Mx is the invariant mass of the diffractive system, Bp is the diffractive slope and W 
the total energy of the 7*p system. 

The dipole cross section, <Tdi P (x,r), is determined by the QCD dynamics, being closely 
related to the solution of the QCD non-linear evolution equations (For recent reviews see, 



e.g. Refs. 3, QQQl) 



(r dip {x,r) 



2 j d 2 bU r {x,r,b) 



(12) 



where is the fundamental dipole-target forward scattering amplitude for a given impact 
parameter 6, which encodes all the information about the hadronic scattering, and thus 
about the non-linear and quantum effects in the hadron wave function. In what follows we 
estimate, for the first time, the diffractive structure function considering the BUW and RC 
BK parameterizations as input in the calculations. Both consider an independence of the 
amplitude with respect to the impact parameter, which implies that o div = a J\f(x,r), <7 
being a free parameter fitted to the respective data. It is important to emphasize that the 
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FIG. 3: Predictions for F 2 D(3) (Q 2 , p, x P ) compared with the ZEUS data Q- Solid line: RC BK 
model; Dashed line: BUW model. 



value of o"o directly constrains the diffractive slope Bp, since both are related by <To = 4tiBd 
if we assume a Gaussian form factor for the proton 291 ] . 

In order to calculate the diffractive structure function and compare with the HERA 
data we need to specify the diffractive slope Bp and the coupling a s , which determine the 
normalization of F 2 D (Q 2 ,/3, xjp). In particular, the magnitude of the qqg contribution is 
strongly dependent on the value of a s . The BUW model implies Bp = 4.3 GeV -2 , while 
the RC BK one implies Bo = 6.7 GeV -2 . Both values are in reasonable agreement with 
the experimental data [81|. On the other hand, we are still free to choose the value of a s . 
Following 65j, we assume a s = 0.15. In a more detailed study we could consider its running 
with Q 2 or perform a fit to experimental data. However, as our goal is to check if these 
models can describe the experimental data, we postpone this study to a future publication. 

In which concerns the RC BK model, it is important to specify the initial conditions for 
the evolution of the scattering amplitude. In 49J, two families of initial conditions were 
considered, one inspired in the GBW saturation model 20j, and the another one in the MV 
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model [56(. In the present analysis, we use the latter, given by 



7V MV (r, Y = 0) = 1 - exp 



L QCD 



4 / VrA 



(13) 



where Q 2 is the initial saturation scale squared and 7 is an anomalous dimension. Both 
parameters are obtained from the fit to F2 data and are given by Q 2 = 0.15 GeV 2 and 
7 = 1.13. 

The diffractive cross section ep — > eXY have been measured by the HI and ZEUS experi- 
ments at HERA tagging the proton in the final state (Y = p) or selecting events with a large 
rapidity gap between the systems X and Y in the case of HI and using the Mx-method 
in case of ZEUS. The distinct methods and experimental cuts used by the HI and ZEUS 
collaborations imply a difference of normalization between their data. Moreover, while the 
ZEUS data are given for the diffractive structure function F 2 D , the HI ones are presented 
for the reduced cross section which is expressed in terms of a combination of diffractive 
structure functions and kinematical factors. As our predictions are for F% , as given in 

821 • Furthermore, as the 



Eq. (El), we restrict our comparison to the recent ZEUS data 
dipole model is more suitable for the description of the diffractive structure functions in 
the region of low and moderate Q 2 , we restrict our comparison to the experimental data 
in the kinematical region of Q 2 < 50 GeV 2 and xp < 10~ 2 . In Fig . [3] we compare the 
predictions of the BUW and RC BK models with the ZEUS data 



for five values of the 



photon virtuality Q 2 . We have that both models describe reasonably the experimental data 
at medium values of /3. However, the predictions differ at large and small f3. In particular, 
the BUW model over predicts the data at large (3 and large Q 2 , while the RC BK model 
describes quite well this same data set. This result can be associated to the features of 
the BUW model, which assumes the geometric scaling property and predicts that the limit 
7 — > 1 is rapidly approached. As a consequence, one obtains different predictions for the 



momentum dependence of the dipole cross section in comparison with, 



or instance, the IIM 



651 ]. It is important to 



model, which is able to describe the data in this kinematical range 
emphasize that in our study the slope parameter is not a free parameter, being derived from 
o"o. An alternative is to assume B^asa free parameter to be obtained from the fit. However, 
this procedure would change the normalization of the BUW predictions by a constant factor 
for all values of (3 and Q 2 . A larger value of Bp would improve the description at large 
P, but would lead to an underprediction of the small (3 data. Consequently, the conclusion 

12 



that the BUW model is only able to describe a restricted kinematical range would not be 
modified. In the case of the RC BK model, it over predicts the diffractive structure function 
at small (3 and Q 2 . Although it is possible to get rid of this over prediction by modifying 
the value of the coupling a s , this would lead to non-physical small values for it. This may 



be an indication that a generalized expression for F% {Q 2 -, (3, Xp), as that proposed in 29]. 
should be used to calculate the diffractive structure function at small-/?. 

IV. FORWARD HADRON PRODUCTION 

Another source of information on QCD dynamics at high parton density is the forward 
hadron production in hadron-hadron collisions at RHIC. In particular, the observed sup- 
pression of the normalized hadron production transverse momentum in dAu collisions as 
compared to pp collisions has been considered an important signature of the Color Glass 
Condensate physics (For a review see e.g. As pointed in Ref. the forward hadron 

production in hadron-hadron collisions is a typical example of a dilute-dense process, which 
is an ideal system to study the small-x components of the target wave function. In this case 
the cross section is expressed as a convolution of the standard parton distributions for the 
dilute projectile, the dipole-hadron scattering amplitude (which includes the high-density 
effects) and the parton fragmentation functions. Basically, assuming this generalized dense- 
dilute factorization, the minimum bias invariant yield for single-inclusive hadron production 



in hadron-hadron processes is described in the CGC formalism by 33|, |3J, [84 1 



d 2 N pp(A)->hX l r I 



dx\- 1 



fq/pfaiPrWr ( x 2, T^PT^j D h/q 



dyd 2 p T (2tt) 2 J Xp x F 

+ fg/ P (xi,PTWA (x 2 , ^Pt^J D h/g C~iPt) 



(14) 



where pt, y and xf are the transverse momentum, rapidity and the Feynman-x of the pro- 
duced hadron, respectively. The variable x\ denotes the momentum fraction of a projectile 
parton, f{x\ : p^) is the projectile parton distribution functions and D{z,p\} the parton 
fragmentation functions into hadrons. These quantities evolve according to the DGLAP 



evolution equations [8g, |8y, |87| and obey the momentum sum- rule. In Eq. ( !T4l) . Afp(x, k) 
and Na(x, k) are the fundamental and adjoint representations of the forward dipole ampli- 
tude in momentum space, which represents the probability for scattering of a quark and a 
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FIG. 4: Inclusive 7r° production cross section in pp collisions at RHIC energies. Data from STAR 
collaboration [88I]. We assume K(r) = 4.0) = K(rj = 3.8) = 1.4 and K(r] = 3.3) = 1.0 for the RC 
BK (solid line) and BUW (long-dashed line) predictions. 

gluon off the nucleus, respectively. Moreover, xp = ^e y and the momentum fraction of the 



target partons is given by £2 = x\e~ 2y (For details see e.g. [841]). 

In the last years, several models have been proposed to describe the hadron spectra in dAu 
collisions, obtaining a satisfactory description of these experimental data 32|, |33|, |34j, |35|, [36] . 
In general, these models have been extended to pp collisions in order to calculate the ra- 
tio Rha without a comparison with the corresponding experimental data. In Ref. ^lj the 
authors estimated the hadron production in pp and dAu collisions in a same theoretical 
formalism and compared these predictions with the experimental data. They fixed the only 
free parameter in the calculations (the .fT-factor) and obtained parameter-free predictions 
for the nuclear modification ratio RhA- The comparison with the experimental data demon- 
strated that the BUW model, which assumes the geometric scaling property, is adequate for 
the RHIC kinematical range. Furthermore, for the first time, it is possible to estimate the 
hadron spectra in pp collisions using as input the solution of the BK evolution equation. In 
what follows we calculate the spectra and compare the BUW predictions with those from 
the RC BK equation. 

For the initial conditions, we use the MV-inspired ones, as it was done in the case of 
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diffraction. In Eq. (ll4p . both fundamental and adjoint scattering amplitudes enters in the 
expression, but only the fundamental one is described by BK equation. Thus, we must 
perform the transformation from Mf to Ma described in the Introduction, but now at the 
level of the x variable, and not the squared saturation scale, since we have been working 
directly with the solution of BK equation. More specifically, from the transformation Q 2 S — > 
Q 2 s Cf/Ca = 4/9 Q 2 , it is straightforward to find the transformation for x through the 
relation Q 2 = (x /x) x . 

In Fig. H] we present our predictions for the production of forward tt° mesons in pp 



collisions at forward rapidities and compare our results with the STAR data 88]. In our 
calculations we use the CTEQ5L parameterization 89j for the parton distribution functions 
and the KKP par ameterization for the fragmentation functions [90]. As in previous calcu- 



lations 1 3c 



34 



35 



36 



4l| there is only one free parameter in our calculation: the A-factor. 



It is determined in such way to provide the best description of the experimental data and 
is fixed for each rapidity. We can see that both models describe quite well the experimental 
data assuming an identical A-factor. This result is directly associated to the similar be- 
havior of the dipole scattering amplitude at large pair separations predicted by these two 
models (See Fig. [1]), which is the region probed in pp collisions at forward rapidities. 



V. CONCLUSIONS 



Recent calculations of the running coupling corrections to the BK equation allowed a 
global fit to the proton structure function to be performed within an unprecedented level 
of precision. This result motivates the study of other observables in order to test the non- 
linear small- x evolution. In this paper we have extended the predictions of the running 
coupling BK equation to diffractive observables in ep collisions and to the forward hadron 
production in pp collisions. Our results demonstrate that the current experimental data 
can be described using the solution of the BK equation. Consequently, we believe that this 
framework is adequate to calculate the observables which would be studied in the future 
colliders. 
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